Primordial Black Hole Formation 
during the QCD Epoch 



Karsten Jedamzik 
Physics Research Program 
Institute for Geophysics and Planetary Physics 
University of California 
Lawrence Livermore National Laboratory 
Livermore, CA 94550 



We consider the formation of horizon-size primordial black holes (PBH's) from pre- 
existing density fluctuations during cosmic phase transitions. It is pointed out that the 
formation of PBH's should be particularly efficient during the QCD epoch due to a sub- 
stantial reduction of pressure forces during adiabatic collapse, or equivalently, a significant 
decrease in the effective speed of sound during the color- confinement transition. Our con- 
siderations imply that for generic initial density perturbation spectra PBH mass functions 
are expected to exhibit a pronounced peak on the QCD-horizon mass scale ~ 1M . This 
mass scale is roughly coincident with the estimated masses for compact objects recently 
observed in our galactic halo by the MACHO collaboration. Black holes formed during the 
QCD epoch may offer an attractive explanation for the origin of halo dark matter evading 
possibly problematic nucleosynthesis and luminosity bounds on baryonic halo dark matter. 
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Ever since the early works by Zeldovich, Novikov [I]], and Hawking it is clear that 
only moderate deviations from a perfectly homogeneous Friedman universe can lead to copi- 
ous production of PBH's at early epochs. Initially super-horizon size overdense regions can 
collapse and convert into PBH's when they enter into the particle horizon. Once a fluctu- 
ation passes into the particle horizon it's subsequent evolution is essentially a competition 
between pressure forces and gravity. When the equation of state is hard, p = p/3 where p 
is pressure and p is energy density, fluctuations with overdensities exceeding a critical value 
(Sp/p) ^ S^ D ~ 1/3 are anticipated to form black holes whereas fluctuations with overden- 
sities less than this critical value are expected to disperse due to pressure forces @ . Here 
overdensities are specified at fluctuation horizon crossing in uniform Hubble constant gauge. 
A typical mass for PBH's formed during a radiation dominated era (e.g. p = p/3) is of the 
order of the horizon mass. The horizon mass is given by 

M H (T) « ImJ — - — VV^£LV 5 , (1) 
HK ' °U00MeV/ V10.75/ ' V ; 

where T is cosmic temperature and g e ff are the effective relativistic degrees of freedom 

contributing to the Hubble expansion. The only numerical, albeit schematic, simulation of 

PBH formation to date 0] indicates that PBH's form with masses somewhat smaller than 

the horizon mass. It has been emphasized that black hole production from pre-existing 

adiabatic fluctuations can be very efficient when there is a period during the evolution of 

the early universe for which the equation of state is soft (p ~ 0) P, |5[]. In the seventies the 

possible overabundant production of PBH's had also been used to rule out a prolonged soft 

Hagedorn equation of state at T ~ lOOMeVp. Nevertheless with the discovery of quarks 

such speculative "dust" -like eras were believed to only occur at temperatures in excess of 

the electroweak breaking scale (T > lOOGeV). 

There are also various schemes for the spontaneous, phase-transition dynamics related, 
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generation of density perturbations on sub-horizon scales during cosmological first order 
transitions and the concomitant production of PBH's. In the context of the QCD confine- 
ment transition Crawford & Schramm J7|] have argued that the long-range color force could 
lead to the generation of sub-horizon density fluctuations which may turn into planetary 
sized black holes even though the details of this mechanism remain to be explored. Hall 
& Hsu || proposed the formation of PBH's during a first order QCD transition from im- 
ploding, supercooled quark-gluon plasma bubbles. Katalinic & Olinto point out that the 
possible tendency of extreme focusing of inwardly traveling sound waves in a quark-gluon 
plasma bubble, reminiscent of the process of sono-luminescence, can lead to density fluc- 
tuations at the center of the bubble sufficiently large for the production of PBH's. Note 
that these schemes predict PBH masses which are far below the QCD-horizon mass. In the 
following we will show that the universe effectively has a soft equation of state during a 
first-order QCD transition which, depending on cosmic initial conditions, may lead to the 
abundant formation of PBH's on the QCD horizon mass scale. 

A first-order color- deconfinement QCD phase transition is characterized by the coex- 
istence of high-energy density quark-gluon phase (quarks and gluons with strong mutual 
many-body interactions) with low-energy density hadron phase (mostly color-singlet pions 
for a gas at vanishing chemical potential) at a coexistence temperature of the order of 
T c ~ lOOMeV [H], jn], Within the simplistic bag model pressure p, energy density p, and 



entropy density s of the hadronic (h) and quark-gluon (qg) phases may be written as |L2], [13 



p h {T) = \ gh ^T* ■ p qg (T) = \g qa ^ - B ■ (2) 



71-2 -A . frjiS . 



p h (T)=g h -T i ; p qg (T) = g qg —T + B ; (3) 

s h (T) = p h ^T 3 ; s qg (T) = \gj^ . (4) 
Here gt ~ 17.25 and g qg ~ 51.25 are the statistical weights of relativistic particles for the 
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individual phases, including the contributions of pions and leptons for the hadronic phase and 
quarks, gluons, and leptons for the quark-gluon phase. The bag constant B acts effectively 
as a vacuum energy density for the quark-gluon phase and accounts for the strong mutual 
interactions of quarks and gluons. For a first-order transition at coexistence temperature T c , 
the conditions of thermodynamic equilibrium are the equality of pressure p = —(dE/dV)s 
and temperature T = (dE/dS)y between hadronic phase and quark-gluon phase. 

One may consider a region of quark/hadron matter sufficiently larged to include material 
in both phases. The pressure, average energy density (p), and average entropy density (s) 
of such a region of quark/hadron matter are 

p = p qg (T c ) =p h (T c ) , (5) 

(p) = fqgPqg( T c) + (1 - fqg)ph(T c ) , (6) 

(*> = /«a«(Tc) + (1 - fMT c ) , (7) 

where f qg is the fraction of space permeated by quark-gluon phase. An adiabatic compres- 
sion of quark/hadron matter in a state of equilibrium phase coexistence induces a conversion 
of low-energy density hadron phase into high-energy density quark-gluon phase (f qg rises). 
During this process the average energy density increases while pressure and temperature 
remain constant. The pressure response of quark/hadron matter to slow adiabatic expan- 
sion, compression, or collapse is therefore negligible. This may be expressed by defining an 
effective speed of sound for quark/hadron matter in a state of phase mixture, such that in 
thermodynamic equilibrium 

holds exactly. Eq. ||] implies that the universe is Jeans unstable during the quark/hadron 
transition for scales much smaller than the horizon length. In contrast, the Jeans length 
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during "ordinary" radiation dominated eras (vs = 1/ v^3) is of the order of the horizon 
length. 

One may wonder if during a cosmic QCD phase transition thermodynamic equilibrium, 
in particular, constant pressure and temperature are maintained. Consider, for example, a 
region of quark/hadron matter at T c . Upon compression or collapse such a region could, in 
principle, superheat which would yield a pressure response, such that v s * 7^ 0. Superheated 
quark/hadron matter may cool by either the growth of existing quark-gluon phase or the 
nucleation of new quark-gluon bubbles. One may estimate the amount of superheating 
77 = (T — T c )/T c at which heating due to adiabatic collapse is balanced by cooling due to the 
nucleation of critically-sized quark-gluon bubbles 



(p) a 3 T? ( 16tt o 3 \ 
^ ~ 30—4 exp — „ — 



where L, and a are Hubble time, latent heat of the transition, and surface free energy 
of the phase boundary, respectively, and we have assumed 77 -C 1. The right-hand-side of 
Eq. [| gives the cooling rate from nucleation of new phase with the exponent the change in 
free energy due to the spontaneous appearence of a critically-sized bubble of quark-gluon 



phase divided by the temperature [13, 14], whereas the left-hand-side of Eq. [] is simply the 
rate of change in energy density during collapse in the absence of cooling. This latter rate 
is approximately given by the energy density of the mildly non-linear fluctuation ~ (p) over 
the gravitational collapse time scale ~ tjj. Quite independent of the prefactors cooling is 
efficient when the exponent in Eq. |^ is approximately 10-20, which yields rj ~ a 3 ! 2 j LT^ 2 . 
Typical parameters of the phase transition have been estimated to be approximately a ~ 
0.02 - 0.1T C 3 H U, 0, U], L w 2 - 15T C 4 [Tj| [21], || with T c w lOOMeV. This implies 



that quark/hadron matter can not sustain superheating by more than 77 ~ 10 3 — 10 5 during 
adiabatic collapse. 
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These considerations illustrate that the pressure response of quark/hadron matter to 
adiabatic compression, or equivalently the effective speed of sound, is dependent on the am- 
plitude and time scale of the compression process. For example, small-amplitude, (5T/T)<r), 
sound waves in pure hadronic phase, or quark-gluon phase, at T c will propagate with the or- 

1 /2 I — 

dinary speed of sound vs = (dp/dp)$ ~ l/v3. However, the pressure response of "mixed" 
quark-gluon and hadron phase to significant adiabatic collapse, (<5(p)/(p)) ~ 1, will be re- 
duced substantially by a factor ~ rj when compared to the pressure response of "ordinary" 
relativistic matter to collapse. Note that the anamoly in the speed of sound during a first- 
order QCD transition, which was independently discovered by Schmid et al. [p4] , may have 
interesting implications for the growth of sub-horizon size initial density perturbations [24 . 

The universe has effectively a soft equation of state during the QCD transition when 
considering the pressure response to gravitational collapse of density fluctuations. It is 
interesting to note that the evolution of average energy density (p) as a function of scale 
factor R during a first-order QCD transition is given by 

(p){R) = (^)\p q9 + \ph\-\p h , (io) 

which yields for the time evolution of R approximately R(t) ~ t 2//3 , akin to a dust-like era. 
Nevertheless, the duration of this dustlike phase is brief. By using the bag model we may 
derive for the ratio of scale factors at the beginning of the transition, Ro, to the value at 
completion of the transition, Ri, from conservation of entropy {R\/Rq) = (g qg /gh)^ ~ 1-44 , 
which implies that the duration of the transition is of order of the Hubble time at the QCD 
epoch. The duration of the transition may be even substantially shorter than this estimate 
if the latent heat is smaller than the bag model value, L « 15T C 4 , as indicated by most lattice 



QCD simulations pL EO, Bfl, 22 



Moderate amplitude, overdense fluctuations entering into the horizon during the QCD 
transition will experience an almost complete reduction in dispersing pressure forces for 
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approximately the duration of the transition. Such fluctuations will have a fraction of a 
Hubble time longer to collapse unhindered by pressure forces than their counterparts passing 
into the horizon during "ordinary" radiation dominated eras. Therefore the demand on 
fluctuation amplitude 5^ CD for successful formation of a black hole during the QCD transition 
should be lessened when compared to the analogous quantity 5f D , 

For almost scale-invariant initial adiabatic perturbations of the Harrison- Zeldovich type, 
which have approximately equal amplitudes at horizon crossing on all scales, PBH's would 
form more abundantly during the QCD epoch than during radiation dominated epochs lead- 
ing to a peak in the PBH mass function at Mbh ~ 1M Q . Note that a peak in PBH mass 
functions on a given mass scale may also result from non-trivial initial fluctuation spectra 



which, however, require fine-tuning of initial conditions [^, ^] . 

For first order transitions which are much longer than the QCD transition the computa- 
tion of an analogous critical overdensity 5 C may be well approximated by using a spherical 
top-hat model for the evolution of the fluctuation , since there are no pressure forces 
between fluctuation and environment. For such transitions 5 C may be substantially below 
£R.d £ or non _ ro tating fluctuations. However, in the case of a first-order QCD transition we 
only expect mild (order unity) decrease in 5® CD when compared to 5^° since the epoch 
of phase coexistence is short. For successful PBH formation during the QCD era fluctu- 
ation and environment may well exist in different regimes during the pre-PBH formation 
evolution, for example, the overdense fluctuation may be in a quark- gluon/hadron mixture 
in phase coexistence at T c whereas the environment may be purely in hadronic phase at 
T<T C . In this case the reduction in pressure gradients is only partial. Nevertheless, since 
pressure gradients are always reduced a decrease in S^ CD compared to 5^° seems inevitable. 
A reliable estimate of £^ CD would have to be obtained with the help of a numerical general- 
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relativistic hydrodynamics code, using an appropriate metric which asymtotically approaches 
the Robert son- Walker metric in regions far away from a fluctuation, and under the inclusion 
of detailed modeling of the QCD equation of state, in regimes at, but also below and above, 
the transition point. 

We have so far restricted our attention to the QCD transition and furthermore implicitly 
assumed that color-confinement in the early universe proceeds via a first order transition. 
Lattice gauge simulations are still not conclusive as to the order of the QCD transition, 
mainly due to finite resolution effects and the difficulties associated with simulating bare 



quark masses |20], [27], |28|. We wish to stress that a possible reduction of the effective speed 
of sound may be a generic feature of cosmic phase transitions and may not necessarily 
be tied to the character of a transition. In fact, a reduction in the speed of sound of 



order 10-20% for a few Hubble times does occur during the cosmic e + e -annihilation 1 29 1 



(the nature of the e + e~-annihilation is very different from that of a first-order QCD color- 
confinement transition). This effect of pressure reduction during phases of e + e~-creation 
and -annihilation is also well known in stellar evolution calculations, commonly referred to 
as the pair-instability, and relates to the conversion of relativistic energy density (photons) 
to e + e~- rest mass energy density. We note that there is the possibility of an enhancement in 
PBH formation on the e + e~-annihilation horizon mass scale of approximately, M ~ 1O 5 M . 

PBH's formed during the QCD epoch (or any other early epoch) may contribute a sig- 
nificant fraction Qbh to the closure density today if only a tiny fraction e of the radiation 
energy density at the QCD epoch is converted into black hole mass density 



Qbh = 5.8 x 10 7 e(T) ( — ( ) h~ 2 . (12) 



.lOOMeV/ U0.75, 

Here h is the Hubble constant in units of 100 km s~ x Mpc -1 . This is because radiation 
energy density redshifts as 1/R A whereas black hole mass density redshifts as 1/R 3 during 
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the subsequent expansion of the universe. Production of PBH's during the QCD epoch would 
also lead to the spontaneous generation of isocurvature perturbations on super-horizon scales, 
even though this isocurvature component is not expected to play a role in the formation of 
large-scale structure unless M BH > 1O 5 M |JU], |3T|. It is, however, interesting to note that 
PBH's form on the peaks of the underlying adiabatic perturbations and PBH number density 
is therefore strongly correlated with the adiabatic density fluctuations. 

The efficiency e(T c ) for PBH formation during the QCD transition can be obtained from 
the statistics of the initial density perturbations and from <5^ CD . It is given by the fraction 
of QCD horizon volumes which are overdense by more than c^ CD , 

POO 

e(T c ) = J sQCD f(S,T c )dS . (13) 

Here f(S, T) is the probability distribution to find a horizon volume overdense by 5 = (5p/p), 
normalized such that J^fdd = 1. Currently favored mechanisms for the generation of 
primordial density perturbations involve quantum fluctuations of scalar fields which drive 
an extended inflationary period of expansion in the very early universe (T^lTeV). Many of 
such models predict a Gaussian probability distribution 

1 = ^^m exp (-\^M)) ■ (14) 

with approximate variance 

Ti/T (l-n) 

^^ Xr '( 5 x,W 6 i.- ) ' ' (15) 

where n is a spectral index [p^ ]. Simple inflationary scenarios predict equal perturbation 
amplitudes for radiation and pressureless matter (CDM). The reader be advised that the 
mass scale M in Eq. [15] denotes the horizon mass in CDM only. For the QCD epoch, 
and assuming a closed universe, this mass scale is M^ DM « 3 x 10~ 8 M Q (T/100MeV)~ 3 /i 2 . 
The simplest inflationary models predict scale-invariance, n = 1 [Q. There are, however, 
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inflationary scenarios which predict the generation of blue spectra |33| (n > 1). Using Eq. [T3| - 
Eq. [15] and assuming 5® CD in the range 5® CD = 0.05 — 0.2, one can infer that for spectral 
index in the range n = 1.67 — 1.79 cosmologically significant PBH formation during the 
QCD epoch (f2pBH ~ 1) m ay occur. This estimate asesses the sensitivity of the required 
spectral index for significant QCD black hole formation on the undetermined 5^ CB . We note 
that in Gaussian models a 5® CT) only somewhat smaller than 5^ D results in PBH formation 
essentially only on the QCD scale due to the steep decrease of / with 5. 

A spectral index as large as n — 1.7 is incompatible with observed cosmic microwave 
background radiation (CMBR) anisotropies 0] and spectral distortions [36], which imply 



72^1.5 — 1.6. Moreover, the required spectral index for significant QCD black hole formation 
may even exceed the above estimate when non-Gaussian, skew-negative features, resulting 
for most inflationary scenarios producing blue spectra, are taken into account |37|]. This 
constraint may be circumvented in two ways. First, it may be that the effective spectral index 
increases with decreasing mass scale. Cosmologically significant PBH formation at the QCD 
scale may then still be compatible with the CMBR limits, since those limits are derived on 
scales much larger than the QCD horizon scale. However, in this case the primordial density 
perturbation spectrum has to be such that constraints derived from PBH formation on mass 
scales smaller than the QCD- horizon mass are not violated ||38|| . These limits may be fairly 



stringent, in particular, when the reheating temperature after an inflationary epoch is high 
(Trh 3> ITeV) and/or the inflationary epoch is followed by a prolonged reheating period. As 
a second possibility it may be that the initial density perturbations are non-Gaussian ||9j f|0| 
and conspire to have a skew-positive tail, such that a fraction of approximately ~ 10~ 8 
horizon volumes are overdense by more than 5^ CB . 

Our findings are particularly interesting in light of the recent observations of the MA- 
CHO collaboration |?1| that a significant fraction of the dark matter in the galactic halo 



may be composed of massive compact objects. Even though the statistics of gravitational 

10 



microlensing events is still poor there is accumulating evidence for a sharp cutoff on the 
distribution of masses of compact halo objects, such that essentially all compact halo ob- 
jects have M > O.1M for a standard halo model. The identification of these objects as red or 
white dwarfs may be problematic. Searches for a dwarf population in our halo with the Wide 
Field Camera on the Hubble Space Telescope yield fairly stringent limits on the fraction of 



mass contributed by dwarfs to the halo |42], |43[ . In the case of an abundant halo white dwarf 
population one may also have to address the problem of overproducing metals. 

In summary, we have argued that PBH formation from pre-existing adiabatic fluctuations 
may be particularly efficient during cosmic phase transitions and periods of particle annihi- 
lation due to a softening of the equation of state. In the case of the QCD color- confinement 
transition we have found that PBH's may form abundantly on the mass scale Mbh ~ 1-^©, 
which is surprisingly close to the inferred masses of compact objects recently discovered in 
our galactic halo by the MACHO collaboration A peak in the mass function of compact 
halo dark matter may, in future, be observationally verified by gravitational microlensing 
experiments. PBH's formed during the QCD epoch provide a natural explanation for such 
a peak, and furthermore evade possibly problematic bounds on baryonic compact halo dark 
matter. 

I am grateful to C. R. Alcock, G. M. Fuller, and A. Olinto for encouragement and valuable 
discussions. This work was performed under the auspices of the US Department of Energy 
by the Lawrence Livermore National Laboratory under contract number W-7405-ENG-48 
and DoE Nuclear Theory grant SF-ENG-48. 
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